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The Schrodinger equation for N particles interacting through effective pair 
potentials is derived from the massless Nelson model with ultraviolet cutoffs. 
We consider a scaling limit where the particles are slow and heavy, but, in 
contrast to earlier work 0, no "weak coupling" is assumed. To this end 
we prove a space-adiabatic theorem without gap condition which gives, in 
particular, control on the rate of convergence in the adiabatic limit. 

1 Introduction 

The physical picture underlying nonrelativistic quantum electrodynamics is that of 
charged particles which interact through the exchange of photons and dissipate energy 
through emission of photons. In situations where the velocities of the particles are small 
compared to the propagation speed of the photons the interaction is given through ef- 
fective, instantaneous pair potentials. If, in addition, also accelerations are small, then 
dissipation through radiation can be neglected in good approximation. Instead of full 
nonrelativistic QED we consider the massless Nelson model. This model describes iV 
spinless particles coupled to a scalar Bose field of zero mass. 

The content of this work is a mathematical derivation of the time-dependent Schro- 
dinger equation for iV particles with Coulombic pair potentials from the massless Nelson 
model with ultraviolet cutoffs. The key mechanism in our derivation is adiabatic decou- 
pling without a spectral gap. 

Before we turn to a more careful discussion of the type of scaling we shall consider, 
notice that the coupling of N noninteracting particles to the radiation field has three 
effects. 



• The effective mass, or more precisely, the effective dispersion relation of the parti- 
cles is modified. The term "effective" refers to the reaction of the particles to weak 
external forces. The physical picture is that each particle now carries a cloud of 
photons with it, which makes it heavier. 

• The particles feel an interaction mediated through the field. If the propagation 
speed of the particles is small compared to the one of the photons, then retarda- 
tion effects should be negligible and the interaction between the particles can be 
described in good approximation by instantaneous pair potentials. 

• Energy is dissipated through photons moving freely to infinity. The motion of the 
particles is, in general, no longer of Hamiltonian type. The rate of energy emitted 
as photons is proportional to the acceleration of a particle squared. 

The scaling to be studied is most conveniently explained on the classical level. The 
classical equations of motion for N particles with positions qj, masses rrij and rigid 
"charge" distributions pj coupled to the scalar field 4>(x, t) with propagation speed c are 



One can think of Pj(x) = Cj ip(x) as a smeared out point charge Cj with a form factor 
if > satisfying J R3 dx <f(x) = 1. Taking the limit c — > oo in ([[]) yields the Poisson 
equation for the field and thus, after elimination of the field, (§) describes N particles 
interacting through smeared Coulomb potentials. Mass renormalization for the particles 
in not visible at leading order. 

Instead of taking c — > oo one can as well explore for which scaling of the particle 
properties one obtains analogous effective equations. Since retardation effects should be 
negligible, the initial velocities of the particles are now assumed to be 0(e) compared 
to the fixed propagation speed c = 1 of the field, £ C 1. In order to see motion of 
the particles over finite distances, we have to follow this dynamics at least over times 
of order 0(e~ l ). To make sure that the velocities are still of order 0(e) after times of 
order 0(e~ l ), the accelerations must be at most of order 0(e 2 ). The last constraint 
also guarantees that the energy dissipated over times of order 0(e~ l ) is at most of order 



The natural procedure would now be to consider such initial data, for which the 
velocities stay of order 0{e) over sufficiently long times. The problem simplifies if we 
assume, as we shall do in this work, that the mass of a particle is of order 0(e~ 2 ). As a 
consequence accelerations are and stay of order 0(e 2 ) uniformly for all initial conditions. 
In this scaling limit mass renormalization is not visible at leading order. Indeed, if we 
substitute t! = et and m'j = e 2 rrij in ([1]) and (g), we find that the limit e — > is equivalent 
to the limit c — > oo. 
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After quantization, however, the two limiting procedures are no longer equivalent. 
The limit c —>■ oo for the Nelson model was analyzed by Davies |5| and later also by 
Hiroshima [§, who removed the ultraviolet cutoff. A comparison of their results with 
ours can be found at the end of this introduction. We will adopt the point of view that 
it is more natural to explore the regime of particle properties which gives rise to effective 
equations than to take the limit c — ^ oo. 

The deeper reason for our choice is that the more natural procedure of restricting 
to appropriate initial conditions gives rise to a similar mathematical structure. If the 
bare mass of the particles of order 0(1), then the proper scaling which yields effective 
equations with renormalized masses was introduced and analyzed for the classical Abra- 
ham model by Kunze and Spohn, see |T(J [19j and references therein. Denoting again the 



ratio of the velocities of the particles and the field as e, they consider charges initially 
separated by distances of order 0(e~ 2 ) in units of their diameter. Hence the forces are 
0(e A ) initially. For times up to order 0(e~ 3 ) and for appropriate initial conditions - 
excluding head on collisions - the separation of the particles remains of order 0{e~ 2 ) 
and thus the velocities remain of order O(e). In particular, the rescaled macroscopic po- 
sition q'(t') = e~ 2 q(t'/e 3 ) satisfies (d/dt') 2 q'(t') = 0(e A ), which matches the order of the 
forces. As a consequence one obtains a sensible limiting dynamics for the macroscopic 
variables. 

One would expect that the same scaling limit applied to the quantum mechani- 
cal model yields in a similar fashion effective dynamics with renormalized dispersion. 
However, inserting this scaling into the massless Nelson model, one faces mathematical 
problems beyond those in the simpler m = 0(e~ 2 ) scaling. Without going into details 
we remark that the main problem is that for massless bosons the Hamiltonian at fixed 
total momentum does not have a ground state in Fock space, cf. |7|, |!|. (As a conse- 
quence it is not even clear how to translate the result in |2l] for a single quantum particle 
coupled to a massive quantized scalar field and subject to weak external forces to the 
massless case.) Nevertheless, the simpler scaling with m = 0(e~ 2 ) provides at least a 
first step in the right direction, since the mechanism of adiabatic decoupling without 
gap will certainly play a crucial role also in a more refined analysis. 

In the remainder of the introduction we briefly present the massless Nelson model, 
explain our main result and compare it to Davies' "weak coupling limit" ||. 

Up to a modified dispersion for the particles, the following model is obtained through 
canonical quantization of the classical system (|IJ) and (0). The state space for iV spinless 
particles is L 2 (R 3iV ) and as Hamiltonian we take 

N 

where c max is the maximally attainable speed of the particles and m their mass, h — 1. 
As explained before, we consider the scaling limit 

e < 1 with c max = 0(e) and m = 0(e~ 2 ) . (4) 

It might seem somewhat artificial to have a relativistic dispersion relation for the parti- 
cles which does not contain the speed of light, but some other maximal speed c max . This 
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is done only for the sake of simple presentation. We could as well consider the quadratic 
dispersion H p = — ^2j=x im^j ^ or ^ e P ar ticles. However, there would be no maximal 
speed and we would be forced to either introduce a cutoff for large momenta or to change 
the topology in (|17]). While both strategies are technically straightforward by using ex- 
actly the same methods as in |2D| in the context of Born-Oppenheimer approximation, 
they would obscure the simple structure of our result. 

We insert the scaling (|4]) into (|3]) and change units such that the particle Hamiltonian 
is now given through 

TV 

H i = E v ~ £2 ^ + 1 • (5) 

3=1 

The particles are coupled to a scalar field whose state is an element of the bosonic Fock 
space over L 2 (R 3 ) given as 

-F=©~ =0 ^)L 2 (M 3 ), (6) 

where is the m-times symmetric tensor product and ®^L 2 (R 3 ) := C. The Hamil- 
tonian for the free bosonic field is 



fr, = dr(|fc|), (7) 

where k is the boson momentum. In our units the propagation speed of the bosons is 
equal to one. The reader who is not familiar with the notation is asked to consult the 
beginning of Section 3, where the model is introduced in full detail. 

In the standard Nelson model the coupling between the j th particle and the field is 
given through 

Hi :j = dy <p(y) pj(y - xj) , (8) 

where <fi is the field operator in position representation and Xj the position of the j th 
particle. The charge density pj G L 1 (IR 3 ) R L 2 (IR 3 ) of the j th particle is assumed to be 
spherically symmetric and its Fourier transform is denoted by pj. For the moment we 
also assume an infrared condition, namely that 

E|i^ 2 (K 3 ). (9) 

3=1 1 1 

Condition (^) constrains the total charge of the system but not that of an individual 
particle to zero. The state of the combined particles + field system is an element of 

H = L 2 (R 3N )®F 

and its time evolution is generated by the Hamiltonian 

N 

h £ = H;®i + i®ar(\k\) + J2 H i,i- ( 10 ) 

3=1 
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Figure 1: The spectrum of Hq(x) for N = 2. The thick line indicates the eigenvalue 
E (x) sitting at the bottom of continuous spectrum. 

Note that H contains no terms which directly couple different particles. All interactions 
between the particles must be mediated through the boson field. 

Our goal is the construction of approximate solutions of the time dependent Schrodin- 
ger equation 

ie—^{t) = H E ^{t) ) $(0) = * £H (11) 

tilL 

from solutions of an effective Schrodinger equation 

ie±,l>(t) = H*gm , #)=^oei 2 (M W ) (12) 

for the particles only. Notice the factor e in front of the time derivative in ( |TT1) and 
(|12"D, which means that we switched to a time scale of order e~ l in microscopic units. As 
explained before, this is necessary in order to see nontrivial dynamics of the particles, 
since their speed is 0{e). 

We remark that the scaling (f|) coincides with the one in time-dependent Born- 
Oppenheimer approximation, where m = 0(e~ 2 ) is the mass of the nuclei and where, at 
fixed kinetic energy, the velocities of the nuclei are also of order 0(e). The Hamiltonian 
( [TOD h as the same structure as the molecular Hamiltonian and the role of the electrons 
in the Born-Oppenheimer approximation is now played by the bosons. 

The key observation for the following is that the interaction Hamiltonian depends 
only on the configuration x of the particles and that the operator 

TV 

H (x) = dT(\k\) + Y,H l , ] (x), 

3=1 

which acts on JF for fixed x G M. 3N , has a unique ground state Q(x) with ground state 
energy 

N j-l 

E o(x) = ^2^2v ij (x i -Xj) + e , (13) 

j=2 i=l 
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where 



V i Az) = - dvdw^— t — 'J^—L 14 



and 

e 



-if / dvdw eM£M. . (is) 

2^VxK3 4vr|t;-w| 

Vij(z) is the electrostatic interaction energy of the charge distributions pi and pj at 
distance z, however, with the "wrong" sign. It is a peculiarity of the scalar field that 
the interaction between charges with equal sign is attractive, eo is the sum of all self 
energies. The remainder of the spectrum is purely absolutely continuous and the ground 
state energy is not isolated, cf. Figure 1. 

Let P*(x) = \Q(x))(fl(x)\, then the states in 

RanP* = {[ dxip(x)n(x) : ip G L 2 (R 3Af )l C H (16) 

correspond to wave packets without free bosons. If the particles are moving at small 
speeds and if the accelerations are also small, one expects that no free bosons are cre- 
ated, i.e. that RanP* is approximately invariant under the dynamics generated by H e . 
Moreover the wave function ip(x) of the particles should approximately be governed by 
the effective Schrodinger equation (O) with 



TV N j-l 

j=l j=2 i~- 



13 3*j j 



Our main result, Theorem [7|, states that for \l/o — J* S dx tpo(x)fl(x) G RanP* we have 

e -iHH/ e ^ o _ j® dx ( e ~iHl s t/e^ ^ e ~ie t/e n ^ = (e^\n{l / s)) (1 + \t\) ||^ || ■ 

(17) 

Notice that in the approximate solution of the full Schrodinger equation the state of 
the field is, up to a fast oscillating global phase e~ ieot ^ e , adiabatically following the 
motion of the particles. In particular, there are no bosons traveling back and forth 
between the particles and the phrase that the particles "interact through the exchange 
of bosons" , which comes from perturbation theory, should not be taken literally in the 
present setting. 

As mentioned before, there is a strong similarity to the time-dependent Born-Oppen- 
heimer approximation, where one obtains an effective Schrodinger equation for the nu- 
clei in a molecule with an effective potential generated by the electrons P0[ . In both 
cases the physical mechanism which leads to the approximate invariance of the subspace 
RanP* is adiabatic decoupling. I.e. the separation of time scales for the motion of the 
different parts of the system lets the fast degrees of freedom, in our case the bosons, 
instantaneously adjust to the motion of the slow degrees of freedom, the particles. 
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However, for massless bosons - in contrast to the Born-Oppenheimer approximation 
- there is no spectral gap which pointwise separates the energy band E$ = {(x, Eq(x)) : 
x G M. 3N } from the remainder of the spectrum of Hq(x), but Eq lies at the bottom 
of continuous spectrum. Hence we need a space-adiabatic theorem, cf. [|(| |2l], |15|] , 



without gap condition. Only recently time-adiabatic theorems without a gap condition 
were established in |2|, |], |22|. The notion time-adiabatic refers to the setting where 
the Hamiltonian of the system is itself time-dependent and varies on a slow time scale. 
In Section 2 a general space-adiabatic theorem without gap condition is formulated 
and proved. The proof is based on ideas developed in and our approach gives, in 
particular, good control on the rate of convergence in the adiabatic limit. 

As an application of the result from Section 2 we consider in Section 3 the scaling 
limit e — > of the massless Nelson model as described above. We emphasize at this 
point that, in view of the missing gap condition, the rate of convergence 0(e^/\n(l/e)) 
in (|TTD is surprisingly fast, since it is almost as good as in the case with a gap. Moreover, 
if all particles have individually total charge equal to zero, then the rate is exactly 0(e) 
as in the case with a gap. Hence, the logarithmic correction must be attributed to the 
Coulombic long range character of the interparticle interaction. In the situation with 
gap it is known [^, [15] that the wave function stays in a subspace RanPf which is 
e-close to the band subspace RanP* up to an error of order 0(e°°). However, in the 
situation without gap, we expect that a piece of order e a , a < oo, of the wave function 
is "leaking out" of RanP* in the sense that it becomes orthogonal to RanP* under time 
evolution. Physical considerations suggest that a = 3/2 for the present problem, see 
Remark P. As a consequence, the e 2 corrections to the effective Hamiltonian are still 



dominating dissipation and can be formally derived from the results in |T5[ . The effective 
Hamiltonian (|55| ) then contains a renormalized mass term and the momentum dependent 
Darwin interaction. 

Finally let us compare our results to those obtained by Davies U , who considers the 
limit c — > oo for the Hamiltonian 

H C = H P ® 1 + 1® dr(c|/e|) + v^Pi- (18) 

Notice that H c is obtained through canonical quantization of (|l|) and (§) if one does not 
set c = 1 as we did before. Davies proves that for all t £ M 

s - lim e- iHH {^ ® Q) = {e~ i{HcS+eo)t ij) <g> Q , 

where Q = {1,0,0, . . .} denotes the Fock vacuum and H p := H p =1 and H e s := iP^f 1 . 
This shows that although the limit c — > oo is equivalent to our scaling on the classical 
level, the results for the quantum model differ qualitatively. While we obtain effective 
dynamics for states which contain a nonzero number of bosons independent of e, cf. 
(|17D, the c — > oo limit yields effective dynamics for states which contain no bosons at 
all. Furthermore, the limit e — > is a singular limit as no limiting dynamics for e = 
exists. 

In Section 3 we also consider an infrared-renormalized model suggested in |l| and 
p~2]1 , which allows us to do without the global infrared condition (|9]). The results are 



exactly the same as in the standard Nelson model. 
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2 A space-adiabatic theorem without gap condition 



Generalizing from the time-adiabatic theorem of quantum mechanics ||, we consider 
perturbations of self-adjoint operators H , which are fibered over the base space R n , 
where, for better readability, we use M := R n to denote this base space. Let Hi be a 
separable Hilbert space and let dx denote Lebesgue measure. Recall that H acting on 



H = L (M, dx) ®7it — L (M, dx; Hi) is called fibered, cf. lfLq| , if there is a measurable 



map M 3 x i— > Hq(x) with values in the self-adjoint operators on Hi such that 

H Q — dx H (x) . 
Jm 

There seems to be no standard name for the set 

S = Ux,s) EM xR,sE a(H {x))} 

and we propose to call it the fibered spectrum of H . 

Let a* C S be such that x i— > P*(x) is measurable, where P*(x) denotes the spectral 
projection of H (x) associated with a*(x). Then P* = J A ® dxP*{x) is an orthogonal 
projection which commutes with H , but which is in general not a spectral projection 
ofP . 

We consider perturbations of Hq, which mix the fibers, in a sense, slowly. As a 
prototype consider for a sufficiently regular real valued function h on "momentum space" 
]R n the self-adjoint operator h £ = h(—ieV x ) on L 2 (M). Here e > is the adiabatic 
parameter and [h £ <g> 1, A] = 0(e) for any operator A which is fibered over M. Let 



H £ = H + h £ <g> 1 



then the invariant subspaces for H constructed above are still "approximately" invari- 
ant for H £ with e small, since [iP,P*] = 0(e) and thus [e" iH£s ,P*] = 0(e\s\). But 
the relevant time scale for the dynamics generated by h £ is t/e with t = 0(1). Thus 
the unitary group of interest is e~ tHH ^ £ . However, according to the naive argument, 
[e~ lHH / £ , P*\ = 0(\t\) and the subspaces RanP,, seem to be not even approximately 
invariant as e — > 0. 

It is well known PUI ETJ that the failure of the naive argument can be cured if a* is 



separated by a gap from the remainder of the fibered spectrum S. Then [e lHH / £ ^ = 
0(e) (1 + \t\), a result that was baptized space-adiabatic theorem in PU| . The object of 



this section is to establish an analogous result without assuming a gap condition. 

We remark that the general setup for space-adiabatic theory are Hamiltonians which 
are "fibered" over phase space, in the sense that they can be written as quantizations of 
operator valued symbols [15| . 



2.1 Assumptions and results 

Let H (x), x G M, be a family of self-adjoint operators on some common dense domain 
V C H{, H{ a separable Hilbert space. Let || • \\h (x) denote the graph norm of H (x) on 



8 



V, i.e., for ip G T>, \\iP\\h (x) = \\Ho(x)ip\\ + H^H- We assume that all the i^o^-norms 
are equivalent in the sense that there is an xq G M and constants C\, C2 < 00 such that 
Ci\\i)\\ Ho (xo) < W\h {x) < C 2 \\ip\\ Ho (xo)- Then 

H = dx H (x) 
Jm 

with domain D(Hq) = L 2 (M) <g> X> is self-adjoint, where here and in the following T> 
resp. D(H ) are understood to be equipped with the || • ||h (xq) resp. || • \\h norm. For 
k G No and £ some Banach space let 

C£(R n ,£) = (/ G C k (W,S) : sup ||fl£7(a:)||g < 00 Va G N n with |a| < fc j . 

£(7ii,7i 2 ) denotes the space of bounded linear operators from Hi to 7i 2 and £ SSL (H) 
denotes the set of bounded self-adjoint operators on 7i. Let | • | be the Euclidean norm 
on W 1 and denote the Hessian of a function A on IR n by V^ 2 'A(x). For the resolvent we 
write R X (A) = (A - A) -1 . Let m > 2. 

Assumption Hg\ Let H (-) G Cg(M, C{V,H t )) and for all x G M let P # (x) be an 
orthogonal projection such that H (x) P*(x) = E(x) P*(x) with P*(-) G C™ +1 (M, C(Hf)) 
and E(-) eCT(M,R). 

In addition one of the following assertions holds: 

(i) For 1 < j < n 

lim esssup || 5 Rm x y iS (H (x)) {d Xj P*){x)P*{x)\\ C (H,) = . (19) 
xeM 



(ii) T/iere is a constant So > and a function 77 : [0, 5q] — > [0, 5o] with r](S) > S and a 
constant C < 00 such that for S G (0, So] and 1 < j < n 

esssup || R E(x) _ i5 (Ho(x)) (d x .P*)(x)P*(x)\\ c{H{) < C S~ l V (S) . (20) 

x&M 



iii) In addition to j\2(\ ) for 1 < k,j < n also 



ess sup 

x&M 



d Xk (R E{xyi s(H (x))(d Xj P*)(x)P*(x)) < CS-'rjiS) (21) 



holds. 



A few remarks concerning Assumption H™ are in order: 

• It is not assumed that P*{x) is the spectral projection of H (x) corresponding to 
the eigenvalue E(x). However, (0) holds pointwise in x whenever P*(x) is the 
spectral projection and has finite rank, cf. Proposition |2|. 
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Inequality (pO| ) is always satisfied with rj(5) = 1. For Assumption H™ (ii) and 
(iii) to have nontrivial consequences on the rate of convergence in the adiabatic 
theorem, rj(5) must satisfy lim^o??^) — 0. These assumptions might look rather 
artificial at first sight, but turn out to be very natural in the proof and also in 
our application. For the simpler time-adiabatic setting, which gives rise to similar 



conditions, we refer the reader to [22 



• The regularity of P*(x) has to be assumed, since it does not follow from the regu- 
larity of H (x) without the gap condition, even if P*(x) is spectral. The regularity 
of E(x) follows from the one of Hq(x) and P*(x) whenever P*(x) has finite rank, 
as can be seen by writing E(x) = tr(if (a ; )P*(a ; ))/trP*(a;). 

The "band subspace" RanP* defined through P* = f® dxP*(x) is invariant under the 
dynamics generated by H , since [Ho,P*\ = holds by construction. We will consider 
perturbations h £ of Hq satisfying 

Assumption h m . For e £ (0, 1] let h £ be a self-adjoint operator with domain D(h) C Ti 
independent of e such that Hq + h £ is essentially self-adjoint on D(h) fl D(Hq). There 
exists an operator (Dh) £ £ £ sa (7^)® n with sup e w ,i] II I (Dh) £ \ \\c(H) < 00 satisfying: 

(i) There is a constant C < oo such that for each A £ C™(M, C(H{)) 

in 

\\[h e ,A]+ieV x A- (DhY\\ C{H) < Cj2 £j SU P 119^0*0 lUcw,) ■ 

j=2 x<EM, \a\=j 

(ii) There is a constant C < oo such that 

|| | [(Dhf,H ) I \\c(D {Ho ),H) + II I l(Dh) £ ,h £ ] | \\ C( H) <eC. 

By assumption, H £ = H + h £ is essentially self-adjoint on D(h) (1D(H ) and we use 
its closure, again denoted by H £ , to define for t £ R 

U £ (t) = e - iHH/£ . 

Since, according to Assumption h m (i), [H £ ,P*\ = [h £ ,P*\ = 0(e), the naive argument 
gives [U e (t),P*\ = \t\0(l). Indeed, our aim is to cure the failure of the naive argument 
and to show that RanP* is invariant for U £ (t) in the limit e —> 0. To this end we will 
compare U £ (t) with the unitary group generated by 

#!iag = H + P*h £ P* + P^h £ P^. 

Also tff iag is self-adjoint on D(H £ ) since P*(-) £ C™(M,£(H f )) and thus H £ - H £ iSLg = 
P^[/i £ ,P*]P* - P*[h e , P m ]Pj- is bounded according to h m (i). A gain we abbreviate for 
t £ R 

U £ diag (t) = e-^WA , 
and we have by construction that 

[A,£/ d £ iag (*)] = 0, 

i.e. RanP* and RanP^ are invariant subspaces for the dynamics generated by H^ iag . 
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Theorem 1. Assume H™ and h m for some m > 2. Let e £ (0, So], then 



Hm 




(i) implies that for t £ 



lim 



U £ (t)-U^ g (t) 



C(H) 



(ii) implies that for some constant C < oo and all t £ 



^(*)-^li ag W < C V {ei){l + \t\) 



C(H) 



Hm 




(iii) implies that for some constant C < oo and all t £ 



^(*) - tfLg(*) _ < C77(£)(l + |t|) 



(22) 



(23) 



(24) 



Note that in Theorem [I| the whole spectrum of possible rates of convergence between 
o(l) and 0(e) as in the case with gap is covered. The estimates for the massless Nelson 
model as an application of Theorem [l] will show that, in principle, all rates can occur. 

The following proposition shows that, assuming the first part of H™ but neither (i), 
(ii) or (iii), then Assumption H™ (i) always holds pointwise in x if P#(x) is the spectral 
projection and has finite rank. The proof is analogous to the one of Lemma 4 in 0. 

Proposition 2. Assume Hq without (i), (ii) or (iii). If P*(x) is the spectral projection 
of Hq(x) corresponding to the eigenvalue E(x) and has finite rank, then 



lim II S R 

8^0 



E(x)-i8 



{H {x)){V x P*){x)P*{x) 



0. 



(25) 



Proof. Since P*(x) has finite rank, the uniform statement ( f25|) follows if we can show 
that lim 5 ^ || S R E (x)~is{H (x)) = for all ip £ Reai(V x P*)(x)P*(x). We have 



lim || i SR E{x) _ iS (H (x)) = lim / p^dX) 
s^o 8^0 j R 



(\-E(x)) 2 + S 2 



f^(E(x)), (26) 



where p,^ denotes the spectral measure of H (x) for ip. Since P*(x) is the spectral 
projection on {E(x)} and since, according to (p9|), Ran(V x P*)(x)P*(x) C RanP t ± (a;) we 
have p 1 p(E(x)) = and thus (f25|). □ 



It is clear from ( PB"|) that additional information on the regularity of the spectral 
measure p^, provides some control on the rate of convergence in (p6|). E.g., if p^(d\) = 
p^(\)d\ with p^ £ L°°(R,d\), then 



p^(d\) 



S 2 



(A - E(x)) 2 + S 2 



< 



d\ 



(A - E(x)) 2 + S 2 



0{S) 



and hence (|20|) would hold pointwise in x with r](S) = S 1 ! 2 . In a sense, the rate 
0(e^\n(l/e)) for the massless Nelson model (|17|) is a consequence of the relevant spec- 
tral measure having a density p(A) ~ A — E(x). 
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We emphasize that (|25|) for all x 6 M does not imply H™ (i), even in the case of 
compact M. This is because for pointwise convergence to imply uniform convergence 
one would need uniform equicontinuity of a sequence of functions. However, in the time- 
adiabatic setting it is indeed sufficient to have (|25| ) for almost all x G /, where / C R is 



the relevant time interval, see 22 



2.2 Proof of Theorem J] 

We start with the standard argument and find that on D(H £ ) 



u £ {t) 



tf e (t) J' d8±{U°(-8)U, 



e 

diag 



--f/ £ (t) / dsU 



(H £ 



Tie 
diag 



where 



Pj- h £ P, 



7" 

w diag 



adj. 



(27) 



P^P^P^PJP, 

In ( p7f) and in the following "± adj." means that the adjoint operator of the first term 
in a sum is added resp. subtracted. Inserting h m (i) into ( p7|) and the result back into 
one obtains 



\u £ {t) 



^diag(^) | 



C(H) 



(2? 



In 



/ rfs f/ £ (-s) (P, x (V.P*) P* ■ (Dhf P* + adj.) U £ diag ( 
Jo 

|) we also used that 

{V x P*){x) = P, ± (x)(V x P*)(x)P,(x) + adj. 



+ 0(e)\t\ 



C(H) 



(29) 



which follows from (V x P*)(x) = (V x P?)(x) = (V x P # )(x)P # (x) + P*(x)(V x P*)(x). 

The nontrivial part in adiabatic theorems is to show that also the remaining term 
on the right hand side of (^) vanishes as e — > 0. Assuming a gap condition, the basic 
idea is to express the integrand, which is 0(1), as the time-derivative of a function that 



is 0(e) plus a remainder that is 0{e) and integrate, cf. |2l], |2(J. The key ingredient in 
this case would be the operator 



F(x) = R E{x) (H (x)) (V 2 P*)(x) P(x) 



(30) 



which is, according to (|29|) , well defined and bounded if the eigenvalue E(x) is separated 
from the rest of the spectrum of Hq(x) by a gap and if P*(x) is spectral. 

The definition (|30| ) is made to give [H ,F] = P^~ (V^P*) P*. However, in absence 
of a gap (RO) is not well defined as an operator on 7ii and, following pi, we shift the 



resolvent into the complex plane and define 

F s (x) = R E{x y iS {H {x)) P^-ix) {V x P*){x) P*{x) 
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One now obtains 

[H (x), Fs(x)} = P^(x) (V x P*)(x) P*(x) + Y 5 (x) (31) 

with 

Y s {x) = -i6R E{x y iS (H {x)) {V x P*)(x)P*{x) . (32) 

Assumptions H™ (i), (ii) and (iii) each imply that lim^o = 0. To see this recall 

that for A(-) G L°°{M,C{H f )) one has 

PIU(H) = ess sup ||A(a;)||£( W ) . 

Note that for better readability we omit the Euclidean norm | . . . | in the notation and 
understand that ||A|| always includes also the Euclidean norm if A is an operator with 
n components. Thus with fl3"2| ) we can make the remainder in fl3~H ) arbitrarily small by 
choosing 5 small enough. However, for the time being we let 5 > but carefully keep 
track of the dependence of all errors on 5. 

By assumption H (-) G C™(M,£(D,Hf)) and P»(-) G C™ +1 (M, C(Ht)), which im- 
plies F s (-) G C{?(M, C{Ht)® n ) and hence, according to h m (i), 

m 

II [h £ ,F s ] \\ c(n) < C £V ^P || 8%F S || £(K) =: Me, 5) . (33) 

Combining (|3l|) and (|33"D we obtain 

= Pt{V x P*)P* +0(\\Y s \\J 1 (e,5)), (34) 

where in (j34|) and in the following 0(a,b,c, . . .) stands for a sum of operators whose 
norm in J~-(7i) is bounded by a constant times a + b + c + . . .. Defining 

Bs = F s -(Dh) e P* -adj., 

one finds with h m (ii) and /2(5) = ||-^s||£(w,d(h )) that 

[# £ ,P 5 ] = [^F^-^rP^+^-^^^P^+^-p/i^^PJ+adj. 

= Pt{V x P*)P*-{Dh) e + adj. + 0(£,||y 5 ||,A(£,5),£/ 2 (5)). (35) 

Now the integrand in ( p8|) can be written as the time-derivative of 

A 5 (s) = -ieU £ (-s) B 5 U £ (s), 

plus a remainder: 

^4,( S ) = U £ (-s)[H £ ,B s ]U e (s) 

= U £ (-s) (Pj-(V x P*)P*-(DhyP* + &d).)U%s) 

+ <D(e,\\Y 6 \\J 1 (e 1 5),ef 2 (5)). (36) 
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Inserting (|36|) into (^) enables us to do integration by parts, 



\U £ (t)-U^Jt)\ 



C(H) 



< 



< 



f da (J;M*)) U £ (s)Ul mg (s) 



C(H) 



+ \t\0[e,\\Y s \\J 1 (e,6),ef 2 (5) 

< \\As(t)\\ C{ H) + \\As(0)\\c(H) 

+ 



C(H) 



J*dsA s (s) f±U s (-s)U^ s (s)) 

+ \t\o(e,\\Y s \\,f 1 (e,6),ef 2 (6)) 
< Ce(2+|t|)||F,|U (H) + \t\0(e,\\Y s \\,f 1 (e : 5),ef 2 (8) 



(37) 



For the last inequality in ( 57T ) we used that ||As(*)iU(W) — C e\\Fs\\c(H) uniformly for 
tGM and that 

^-U*(t ,s)U° iag (s,t ) = - l -U%t ,s) (H%s) - H^ s (s))U^ g (s,t ) 

is bounded uniformly, according to fl27|) and h m (i). 

Writing out the various terms in (|3"7|) explicitly, we conclude that there is a constant 
C < oo such that 



F £ (0-^ia g (*)|L (w) <C £ ||F 5 || £(H) + C|t| (e+\\Y 5 \\ c{n) 

m 

+e \\Fs\\c { h,d(Ho)) + £ \\F s \\c(H) + £J SU P U3Mk*)) • 



l«|=j 



(38) 



Hence we are left to establish bounds on F$, on its derivatives and on Y$ in terms of 5, 
which is the content of the following Lemma. 

Lemma 3. Assume H™, then F$(-) e C™(M, C(TLi)® n ) and there is a constant C < oo 
such that for 5 G (0, S ] 

\\F s \\c { H,D { Ho)) < jV(S), (39) 



C 



sup || d% F 5 \\ c{n) < — - r)(S) forl<j<m. 

\a\=j 



(40) 



In case H™ (i) holds, we have ( |3^ ) and ( $Q ) with r](5) = 1. Furthermore, if H™ (i) 
/ioWs, t/ien lim^-^o ||^IU(w) = and if H™ (ii) or (iii) holds, then \\Ys\\c(h) < Ci](5). 
//H™ (iii) holds, then (j^J\) can be improved to 



sup || F 5 || £(H) < g V (S) forl<j<m. 



(41) 
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Before we turn to the proof of Lemma |3| we finish the proof of Theorem [I]. Assuming 
Hq 1 (i), ( |2"2"| ) follows by inserting the bounds from Lemma |3] into (38) and choosing 
5 = 5(e) such that \im £ _> 5(e) = and \im £ ^ e/5(e) 2 = 0. 

If Hq 1 (ii) holds, then the bounds fl39f) and fl40|) inserted into (^) yield 



rw-cw*)!^ <o e 2©+o( e+ ,(«) +e na (42) 

In ( fl"2"l) the optimal choice is 5(e) = ea, which gives (p3[) . Finally, the bounds (^) and 
( f4"T|) inserted into (^Sj) yield 

r (*) - ^(OIU < *M + C (e + V (5) +e M + £e>*M) \t\ , 



3=1 

where the optimal choice 5(e) = e gives (|4]). 



Proof of Lemma [|. We abbreviate REf x )-is(H Q (x)) as R(5,x) in this proof and note 
that P(V) e Cg(M : £(Hi)) and thus G C^(M, £(W f ) en ) follow from #<,(•) e 

C^{M,C{V,H f )) together with P*(-) G C™ +1 (M, £(ft f )). 

We start with the case H™ (ii), where (i) is included by making the obvious changes 
for rj(5) = 1. Assumption H™ (ii) immediately yields 

\\F s \\c(H)< jV(S) (43) 

and the bound on Yg. ( |39| ) follows from H (x)R(5, x) = 1 + (P(x) — i5)R(5, x) and ([43|) 
together with the assumption that E(x) is uniformly bounded. 
For ( pE0| ) we start by observing that 

V x P 5 (x) = P(5,x)V*((V x P*)(x)P*(x)) + (y x R(5,x)yV x P*)(x)P*(x) (44) 
= R(5,x)V x ((V x P*)(x)P*(x)) - R(5,x)(V x H (x) - V x E(x))F 5 (x). 



Using ( ff3D and the fact that |V a; £'(a;)| is uniformly bounded by assumption, we infer 
from (|D that 

II \V X F S \ \\c(n) = sup || \V x F s (x)\ \\ C (H f ) < C(5- 1 + 5~ 2 V (5)) . 

Hence (|0]) follows for j = 1, since 5 < rj(5) by assumption. 

By differentiating (53) again, we find, using a reduced notation with obvious meaning, 
that 

V (2) P 5 = -2P 5 (VP - VP) VP* + R 5 V (2) ((VP,)P,) - P 5 (V (2) P - V (2) P) P 5 . (45) 

Hence ||V^P 5 |U( K ) < C(5~ 3 r](5) + 5' 1 + 5- 2 t](5)) which proves © for j = 2. By 
repeated differentiation one finds inductively ( f4"U| ) for j < m. 

To show (|4T|) assuming H™ (iii), note that (fill) holds by assumption for j = 1 and 
inserted into (g5|) it gives || V (2) P 5 ||,c( W ) < C(6~ 2i q(6) + + tf" 2 ^)). Analogously the 
estimates for all larger j < m are improved by a factor of 5. □ 
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3 Effective dynamics for the massless Nelson model 



As explained in the introduction, we consider iV spinless particles coupled to a scalar, 
massless, Bose field with an ultraviolet regularization in the interaction. This class of 
models is nowadays called Nelson's model Jl], [TJJ after E. Nelson [ 14[] , who studied 



the ultraviolet problem. We briefly complete the introduction of the model and collect 
some basic, well known facts. 

A point in the configuration space M. 3N of the particles is denoted by x = (x±, . . . ,xn) 
and the Hamiltonian for the particles is defined in (|5|). is self-adjoint on the 
domain if 1 (M 3iV ), the first Sobolev space. 

The Hilbert space for the scalar field is the bosonic Fock space over L 2 (R 3 ) defined 
in (D). On D(\fN), M the number operator, the annihilation operator a(f) acts for 
/ G L 2 (M 3 ) as 

K/)^)^!, ...,k m ) = V^+T [ dk f(k) ij (m+1 \k, h,..., k m ) , 



where V = (ip {0 \ ip {2 \ . . .) G D(^N) if and only if E™=o m ll^ (m) P < 00 ■ The 
adjoint a*(/), which is also defined on D(vC\/°), is the creation operator and for f,g G 
L 2 (M 3 ) the operators a(f) and a*(g) obey the canonical commutation relations (CCRs) 



W), «*(<?)] = / dk f(k) g(k) =: (/, g) , [a(f), a(g)} = [a*(f), a*(g)\ = . (46) 

It is common to write a(f) = J R3 dk f(k) a(k). The Hamiltonian of the field as defined 
in (^) can formally be written as 



H { = / dk\k\a*(k)a(k). 

More explicitly, on the m-particle sector the action of Hf is (ilf?/>) (m ' ) (/ci, . . . , k m ) = 
Sj=i l^jl^^^i' • • • > k m ), and H$ is self-adjoint on its maximal domain. For / G L 2 (M 3 ) 
the Segal field operator 

*(/) = ^ (<*(/) + <*•(/)) 

is essentially self-adjoint on D{\fM). The field operator as used in (§) is related to 
$ through (f)(f) = ^{f / \f\k\). For the following it turns out to be more convenient to 
write the interaction Hamiltonian in terms of $, where 



H Y = $(\k\v(x,kj\ 



acts on the Hilbert space 7i = L 2 (M. 3N ) <g> T of the full system. We will consider two 
different choices for v(x,k) in more detail. For the standard Nelson model (SN), as 
discussed in the introduction, one has 



W^) = £e^§§. (47) 
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For the infrared- renormalized models (IR), as considered by Arai Jl] and, more generally, 
by Lorinczi, Minlos and Spohn |T2]], one has 



M^) = E(^-i)jp- ( 48 ) 

In both cases, the charge distribution pj G L 1 (M 3 ) of the j th particle is assumed to be 
real- valued and spherically symmetric. As to be discussed below, cf. Remarks [| and || 
we have to assume the infrared condition @ for the (SN) model, but not for the (IR) 
model. The infrared condition implies, in particular, that the total charge of the system 
of N particles must be zero. 

The full Hamiltonian is given as the sum 

H £ = H; ® 1 + 1 ® H { + Hi + V m ® 1 (49) 

and is essentially self-adjoint on D(H p ®l)nD(l®Hf) if sup x || \k\ s v (x, k) \\l 2 (r 3 ) < oo for 
s G {g, 1}- Only in the (IR) model a potential Vir is added, which acts as multiplication 
with the bounded, real-valued function 

M*) = E / dk ^ ( |nf ] V a -^ + ~ [ dk 

Ji L — 1 



\k\ 



(50) 




Remark 4. If the charge distributions satisfy the infrared condition @, then the (SN) 
Hamiltonian and the (IR) Hamiltonian are related by the unitary transformation 



Uq = exp — i$ 



cf. Jl], which is related to the Gross transformation |TJJ for x — 0. If the infrared 
condition is not satisfied, the (SN) model and the (IR) model carry two inequivalent 
representations of the CCRs for the field operators. Physically speaking, the transfor- 
mation Uq removes the mean field that the N charges would generate, if all of them 
would be moved to the origin. The vacuum in the (IR) representation corresponds to this 
removed mean field in the original representation, a fact which has to be taken care of 
in the interaction: for each particle the interaction term is now evaluated relative to the 
interaction at x = 0, cf. fl4"8|) , which makes also necessary the counter terms Vir. If the 
total charge of the system is different from zero, then the mean field is long range and, as 
a consequence, the corresponding transformation is no longer unitarily implementable. 
Indeed, it was shown that the (SN) Hamiltonian with confining potential does not have 



a ground state, cf. [11], while the (IR) Hamiltonian with the same confining potential 



does have a ground state, cf. CI]. <0 
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In order to apply Theorem |] we observe that H\(x) acts for fixed x G ~R 3N (= M) on 
T (= Hi) and with if „ 0*0 = H t + HAx) + tWa?) we have 



H e = H £ p (g)l+ dxH (x) [^h £ + H 

J M 



The following proposition collects some results about H (x) and its ground state. Its 
proof is postponed to after the presentation of the the main theorem. 

Proposition 5. Assume that v(x, ■) G L 2 (M 3 ) for all x G M. 3N and that for some n > 1 

(i) | • |<9>(x, •) G L 2 (M 3 ) /or a// x G K 3Af and < \a\ < n, 

(ii) sup^siv || a/H <9^(x, •) ||l2 (m3) < oo /or < |a| < n, 

(iii) sup xeIR 3iv || d£v(x, •) ||l2 (r3) < oo for 1 < |a| < n. 

Lei lm(v(x, •), V x w(x, O)^^) = for all x G R 3N and V m {-) G C 3 (M 37V ). Then 

1. H (x) is self adjoint onV = D(H { ) for all x G R 3N and #„(•) G C£(R 3Af , C(V,F)), 
where T> is equipped with the graph-norm of H{. 

2. H (x) has a unique ground state Q(x) for all x G M 3N and, in particular, 
H (x)n(x) = E(x)tt(x) for 



E(x) = -\f 



dk\k\ \v(x,k)\ 2 + V 1R (x). (51) 



Furthermore ft(-) G C£(M 37V , T). 

It is straightforward to check that Im(y(a;, •), V x v(x, -))l 2 (r 3 ) = for vsn defined in 
( pE?D and vir defined in (fJSj). For the (SN) model as well as for the (IR) model ( [51] ) is 



easily evaluated and one finds P(x) = E (x) as given in (p~3|) . 

Remark 6. For the (SN) model the assumptions made on v(x,k) in Proposition |5| 
are satisfied, if Pj(k) decays sufficiently fast for large \k\ and each j = 1,...,N, or, 
equivalently, if Pj{x) is sufficiently smooth. This is an ultraviolet condition individually 
for each particle. But v(x, •) G L 2 (IR 3 ) follows from v(0, •) G L 2 (R 3 ), which is exactly 
the global infrared condition (g). 

While the necessity for an ultraviolet regularization remains in the (IR) model, the 
infrared condition is replaced by J^. \k\~^pj{k){e lk ' X: > — 1) G L 2 (R 3 ), which can be satis- 
fied without having X^ Pj(O) = 0- Thus the (IR) model allows us to consider particles 
with total charge different from zero. <0 

Let P t (x) = \Q{x))(Q{x)\, then P*(-) G C^(M 3iV , £(.F)) and RanP* is a candidate for 
an adiabatically decoupled subspace. Indeed, we will show that h £ = H £ (g) 1 satisfies 
Assumption h m with (Dh) e x . = —ieV x J ^—e 2 A X] + 1 and that H and P* satisfy As- 
sumption Hq 1 (iii) with r;(5) = 5yTri(l/5) if particles with charges different from zero 
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are present and rj(S) = 5 if all particles have total charge zero. Hence we can apply 
Theorem PI to conclude that for some constant C < oo 



e -iHH/e _ e -iH^t/e 



< C V (e)(l + \t\), 



(52) 



with tff iag = P* H* P* + (1 - P*) /F (1 - P*). 

Next observe that the ground state band subspace RanP* is unitarily equivalent to 
the Hilbert space L 2 (M. 3N ) of the N particles in a natural way. Let 



U : RanP* — > L 2 (M. 3N ) , ip ^ (Uip)(x) 
then it is easily checked that for (p G L 2 (M 3Ar ) 



(fi(x) s V(x)) 



(53) 



J R 3iv 



Hence the part of P| iag acting on RanP* is unitarily equivalent to 



UP*H £ PM* 



acting on the Hilbert space of the N particles only. The following theorem shows that 
ip; ff has, at leading order, exactly the form expected from the heuristic "Peierls substi- 
tution" argument. 

Theorem 7. Let H e be defined as in ( f^ ) with v(x,k) either vsT<s{x,k) as in or 
viR_(x,k) as in /^j). For 1 < j < N let pj G L 1 (M 3 ) such that pj satisfies \k\ s pj(k) G 
L 2 (M 3 ) for s G { — 1,4}. For the (SN) model assume, in addition, the infrared condition 
§). Let 

N j-1 
j=2 i=l 



with Vij(z) and eo as in ftL$ ) and fclh\). Then there is a constant C < oo such that 



(54) 



where r](e) = e^/\n(l/s). If all charges satisfy the infrared condition individually, i.e. if 
Pj( k )/\ k \ 3/2 e L 2 (R 3 ) for allj = 1, . . . , N, then $5%) holds with 77(e) = e. 

For sake of better readability the global energy shift eo was, as opposed to (|T7|). 
absorbed into the definition of iP; ff . 

Remark 8. The unitary U intertwines the position operator x® 1 on L 2 (R 3N ) <g> T with 
the position operator x on L 2 (R 3N ) exactly and the momentum operator —ieV x ® 1 on 
L 2 {R 3N ) ® T with — ieV x on L 2 (IR 37V ) up to an error of order e. Thus one can directly 
read off the position distribution and approximately also the momentum distribution of 
the particles from the solution ijj(t) = e'^^^ipo of the effective dynamics. It is not 
necessary to transform back to the full Hilbert space using U*. (} 
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Remark 9. From the discussion of the introduction one expects, on physical grounds, 
that the energy lost through radiation is of order 0(e 3 ) after times of order 0(e~ 1 ). 
Hence the error of order 0{eyj\a.l/e) in fl54|) is not optimal in the sense that the er- 
ror does not correspond to emission of free bosons and thus to dissipation. Indeed, we 



expect that the situation is similar to adiabatic perturbation theory with gap, cf. [pL5 
There should be a subspace RanP^ which is e-close to RanP* and for which the analo- 

3 

gous expression to (B^) holds with an error of order 0{e^ ), possibly with a logarithmic 
correction. The corresponding effective Hamiltonian would then contain two additional 
terms of order e 2 , which, for the case of quadratic dispersion 

N 2 

H £ = -Y-A x . 



for the particles, can be calculated using formula (47) in [TjJ . As result we obtain for the 



Weyl symbol of the effective Hamiltonian including the momentum dependent Darwin 
term 

^(P.!) = E/7l? + £ («)4E / dk (55) 

i= l Zm i 1 3 <t J ^ 3 

with rrij — 1/(1 + and 



the electromagnetic mass. As explained above, for the rigorous justification of fl55|) a 
space-adiabatic theorem without gap but for rotated subspaces is needed and thus 
it is beyond the scope of the present paper. (} 

Before proving Theorem |7| we make up for the 

Proof of Proposition^. A standard estimate (cf. e.g. 0] Proposition 1.3.8) shows that 
for / G L 2 (R 3 ) and any a > 



II// 



i-i 



mm < 4m\% + \ — v ' 1 llL2(R3) + 2 )w\%. m 



Hence $(/) is infinitesimally -Pf-bounded whenever ||/||i,2( K 3) + ||//y| • |||l 2 (r 3 ) < oo. 
Then Kato-Rellich implies that H (x) is self-adjoint on P(Pf), since y/\ ■ \v(x, •) G L 2 
is assumed. Using (i), (ii) and Vir(-) G Cb(IR 3Ar ), we obtain from flSED that 

d%H (x) =${\k\ d%v(x, k)) + d a x V lK {x) (57) 

is relatively bounded with respect to Pf for |a| < n. Moreover, (ii), (|56|) and (|57D imply 
that P (-) G C£(R 37V ,£(P,.P)). 
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To compute the ground state energy E(x) observe that from "completing the square" 
one finds 

H (x)= [ dk\k\ (a*(k) + V *^ k L. X>> ) (a(k) + ^S.) - 1 f dk \k\ \v(x, k)\ 2 + V m (x) . 

It is well known that the map a(k) i— > a(fc) + t>(x, k)/y2 comes from the unitary trans- 
formation £/(x) = exp yi$(iv(x, -))J, i.e. 

f/(x) a(/) CT(x) = a(/) + (/, v(x, •)) , (58) 

whenever v(x, •) G L 2 . Equation fl58|) follows from the fact that [A, [A, £?]] = implies 
[exp (iA) , B] = exp(iA)[iA, B] and the CCRs. Transformations of the form ([58]) are 
called Bogoliubov transformations. 

Therefore H (x) = U(x)H f U*(x) + E(x) with E(x) = -\ J* K3 rfA; |fc| |v(x, A;)| 2 + 
Vir(x). Since i/f f2 = for the unique ground state Q = (1, 0, 0, . . .) G -T 7 , we find 

H (x)Q(x) = E(x)Q(x) with Q(x) = U(x)Q ■ 

Next we need to take derivatives of U(x) with respect to x. It follows from the CCRs 
© that for f,ge L 2 (R 3 ) 

W)M9)]=iMf,9), 

and thus, by assumption, that V x Q(iv(-,x)) = $(iV x v(-,x)) commutes with $(iv(-,x)). 
Hence we obtain that on D(\fM) 



V x U(x) = U(x) i$(iV x v(-, x)) = i<3>(iV x v(-, x)) U(x) . (59) 

By further differentiating (|59| ) we can get up to n th derivatives since d^v(x, k) G L 2 (M 3 ) 
for \a\ < n. In particular we find with (iii) that Vi{x) = U(x)fl G C£ (K 3N , T). □ 

Proof of Theorem 0. We start by showing that the assumptions of Theorem |l| are indeed 
satisfied and thus (|52| ) follows. 

It is straightforward to check that the assumptions on pj imply the assumptions of 
Proposition [5] for n = 5. Hence the first part of Hq follows with P*(x) = \Q(x)){Q(x)\. 
For Hq (iii) observe that, using (|59|), V x Q(x) = V x U(x)Q = U(x) i<&{iS/ x v (•, x))Q , 
and thus 

(tt{x), V x tt{x))r = (O , z$(iV x v(-, x))n )r = . (60) 
As a consequence, (V x P*)(x)P*(x) = \V x ft(x))(ft(x)\. Hence we obtain for 1 < j < N 

R{6, x)(V x .PA(x)P*(x) = i \U(x) Rf(6) *(iV XJ v(; x))Q )(Q(x)\ , (61) 

where Rf(5) = (Hf — iS)^ 1 . For (|20|) one therefore finds 

\\R(8,x)(V Xj P*)(x)P*(x)\\ 2 c{r) = \\R { (5)^\/ Xj v(x,-))n \\% 

Qkl-iS^ke^pjik^kl-l 2 ^ (62) 
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Whenever pj satisfies the infrared condition pj(k)\k\~^ G L 2 (IR 3 ), ( |62"D is bounded uni- 
formly in 5 since | id) | < 1. In general we only assume that pj(k)\k\~ l G L 2 (IR 3 ). 
Using that f>j{k) is bounded uniformly according to Riemann-Lebesgue, fl62| ) becomes 

{Ikl-td^ke^pjik^kl"! 2 <d / j^ + C 2 < Cm 1/5, (63) 

for 5 G (0, -]. To obtain the estimate (|2~ll), we differentiate ( |6lD and find 



V^P(5,x)(V. 3 P0(x)P*(x)J =z|?7(z)i? f ( < y)$(iV a!J i;(x,.))no)(V !Bi n(a:)| 
- \U(x) $(iV Si v(x, •)) i?f(5)$ (iV^ufo OWM^OI 
+ z it: f (5)$(iVg^(x, .))fto)<ft(aO| • (64) 

All three terms in (|64D can be bounded by using the same type of arguments as in fl62]) 
and fl6"3"|): For the first term use in addition that || V Xi fi(x)||-^ f < C and for the second 
term one has to estimate the components in the 0-boson sector and in the 2-boson sector 
separately In summary we showed that H 4 (iii) is satisfied with r] (5) = 5^\n(l/5), and 
with 7](5) =5 if all charges satisfy the infrared condition individually. 

We are left to check for Assumption h 4 . Let h £ = H £ <S> 1 = Ylj h(—i£^7 X j) ® 1 3 
with h(p) = \J p 2 + 1. Essential self-adjointness of h £ + Ho on D(H P g) 1) n D(Hq) is a 
standard result, cf. Proposition 2.1 in [|lj. We define (Dh) £ = (Vh)(—ie'V Xj ) <g> 1, with 
|| (D/i) £ || £ ^- ) ©3iv < 1, and postpone the technical proof of the following Lemma to the 
end of this section. 

Lemma 10. h £ and (Dh) £ satisfy h 4 (i) and (ii). 

We conclude that all assumptions of Theorem [l] are satisfied for the (SN) and the 
(IR) model and thus (|52| ) holds. However, ([54]) follows from ( |52"D by the following Lemma 
and an argument like (p7|). □ 



Lemma 11. There is a constant C < oo such that for e > sufficiently small 

\\{H £ diag -WH £ s U) P*\\ c{n) < e 2 C. 

Proof. In order to apply h 4 (i) to U(x), we have to extend U(x) : RanP*(x) — > C defined 
in (H) to a map U{-) G C 4 (M 3Af , £(.F)) first. To this end let W(ae) = |n )(fi(a;)| and 
note that U*U = P*. With this definition one finds 

F| iag P* = H P* + P*h £ P* = EP* + P*U*h £ UP* + P4h £ ,U*}UP* 
= U*H* e UP* + P*[h £ ,U*]UP*, 

and we are left to show that ||P* [h £ , W]UPA\ = 0{e 2 ). Using h 4 (i) with A = W we 
find that 

P*\h £ ,W\UP, = -ieP*(y x W)-(Dh) £ UP* + 0(e 2 ). 
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However, according to (|60|) 

p*(x)(y x u*)(x) = \n(x))(n(x), v x n(x))(n \ = o, 

and thus the desired result follows. □ 

Proof of Lemma [7Z]. Heuristically h 4 (i) and (ii) hold, because they are just special cases 
of the expansion of a commutator of pseudodifferential operators. However, since h e is 
unbounded and A is only 4-times different iable, we need to check the estimates "by 
hand". 

For notational simplicity we restrict ourselves to the case N = 1, from which the 



general case follows immediately. Let g{p) = 1/a/p 2 + 1, g e = g(—ieV x ) (g> 1 and 
A e C 4 (M 3 ,ft f ), then | ■ \ s g e L^M 3 ) for s G {0,4} and thus for tjj e S 

= J dyg{y) ^A(x) - ey ■ V A(x) + e 2 J ds (y,V^A(x — sey)y)j ip(x - ey) 
= (Ag e ^(x) - ie(vA-Vg £ 4>^)(x) 

+ J dyg{y)e 2 J ds (y, V (2) ' A(x - sey) y)ip(x - ey) . (65) 

From (^) one concludes after a lengthy but straightforward computation involving 
several integrations by parts that 

e 2 A x [<f , A] = -i eVA ■ (Vg) £ (e 2 A x ) + R 

with 

4 

pi ii <cyy sup \\d:A(x)\\ cm . 

~2 xeM? N , \a\=j 

Hence we find 

[h £ , A] = [(1 - e 2 A x )g £ , A] = (1 - e 2 A x )[g £ , A] - [e 2 A x , A]g £ = -xeVA ■ (Dhf + R' 
with 

4 

ii#h< cyy sup n^(x)ii £(Wf) . 

j=2 xm 3N , \a\=j 

This proves h (i). By the same type of arguments one shows also h 4 (ii). □ 
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